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Abstract 

The distribution of the 'k^'k~^'k~ decays in the Dalitz plot has been 

measured by the NA48/2 experiment at the CERN SPS with a sample of 4.71 x 10^ 
fully reconstructed events. With the standard Particle Data Group parameterization 
the following values of the slope parameters were obtained: g = (— 21.134ib0.017)%, 
h = (1.848 ± 0.040)%, k = (-0.463 ± 0.014)%. The quality and statistical accuracy 
of the data have allowed an improvement in precision by more than an order of 
magnitude, and are such as to warrant a more elaborate theoretical treatment, 
including pion-pion rescattering, which is in preparation. 
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Introduction 



The K"^ — s> -K^-K^Ti decay can be described [T] in terms of the Lorentz invariant kinematic 
variables u and v defined as 

u = — , v = — , 1 

ml 

Sij = {Pi + PjY, z,j = 1,2,3, i < j] So = ^{si2 + si3 + S23) ^^mj^ + ml. 

Here rriT^ is the charged pion mass, rriK is the charged kaon mass. Pi are the pion four- 
momenta, the indices i,j = 1, 2 correspond to the two pions of the same electric charge 
("even" pions), and the index i,j = 3 to the other ("odd") pion. The experimental 
Dalitz plot distribution of the -—>■ tt^tt+tt^ decay has been up to now experimentally 
analyzed [2], [3] in terms of a polynomial expansion in powers of u and v. 

\M{u,v)\^ ~ C{u,v) ■ (1 + gu + hu^ + kv^). (2) 

Here g, h, k are the linear and quadratic slope parameters (terms proportional to odd 
powers of v are forbidden by Bose symmetry, and the variable v is defined only up to a 
sign), and C{u,v) is the Coulomb factor [4j: 

C{u,v)= Yl {'^ii/le"'' - 1)}, Uij = 2Tiaqiqj/Pij, (3) 

ij=l,2,3; i<j 

where qi = ±1 are the pion charges, a is the fine structure constant, and jSij is the relative 
velocity of the pions i and j, expressed via the squared invariant mass of the pion pair Sij 
as 

Among the measurements of the slope parameters performed in the past, the most 
precise are reported in [2] based on 1.5 x 10^ decays, and in ^ based on 0.225 x 10^ 
decays. 

The primary goal of the NA48/2 experiment at CERN SPS is the search for direct 
CP-violating charge asymmetries of Dalitz plot linear slopes in Tr^Tr+vr" [5] and 

_^ ^±^0^0 ^ decays. However, the large data sample collected also allows a study of 
the Dalitz plot distributions to be performed at a new level of precision, estimating the 
detector acceptance with a detailed Monte Carlo (MC) simulation. 

The parameterization ([2]) of the tt^tt+tt^ decay distribution takes into account 

electromagnetic interactions only in the first approximation [4] and totally neglects pion 
rescattering effects [7j-[.10j, which were recently shown by NA48/2 to contribute signif- 
icantly to the distribution of — > n^n^n^ decays yjj. The current study aims to 
measure the Dalitz plot slopes within the conventional framework ([2]), and represents a 
first step towards a more complete description of the decay distribution. In particular, the 
presence of radiative effects and strong final state interactions, both ignored in the current 
work, imply that the slope parameters obtained, while fully comparable to those defined 
in the Particle Data Group compilation py, should not be attributed any precise physical 
significance as far as the parameterization of the weak decay in terms of more fundamental 
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Figure 1: Schematic lateral view of the NA48/2 beam line (TAX17,18: motorized beam 
dump/collimators used to select the momentum of the and K~ beams; FDFD/DFDF: 
focusing set of quadrupoles, KABESl-3: kaon beam spectrometer stations), decay volume 
and detector (DCHl-4: drift chambers, HOD: hodoscope, LKr: EM calorimeter, HAG: 
hadron calorimeter, MUV: muon veto). Note that the vertical scales are different in the 
two parts of the figure. 



parameters is concerned; this would require an improved theoretical framework which is 
just being developed [8], and its implementation is postponed for a forthcoming analysis. 

NA48/2 collected data during two runs in 2003 and 2004, with about 50 days of 
efficient data-taking in each run. This analysis is based on 4.71 x 10^ fully reconstructed 
K"^ TT^TT+TT^ events, corresponding to 55% of the data sample collected during the 
NA48/2 run of 2003. 

1 Beams and detectors 

A beam line providing two simultaneous charged beams of opposite signs overlapping in 
space all along the decay volume was designed and built in the high intensity hall (ECN3) 
at the CERN SPS. The beam line is a key element of the experiment, as it allows decays of 

and K~ to be recorded at the same time, and therefore leads to cancellation of several 
systematic uncertainties for the charge asymmetry measurement. Regular alternation of 
magnetic fields in all the beam line elements was adopted. The layout of the beams 
and detectors is shown schematically in Fig. [H The setup is described in a right-handed 
orthogonal coordinate system with the z axis directed downstream along the beam, and 
the y axis directed vertically up. 

The beams are produced by 400 GeV/c protons impinging at zero incidence angle on a 
beryllium target of 40 cm length and 2 mm diameter. Charged particles with momentum 
(60 ± 3) GeV/c are selected in a charge-symmetric way by an achromatic system of four 
dipole magnets with zero total defiection ('achromat'), which splits the two beams in the 
vertical plane and then recombines them on a common axis. Then the beams pass through 
a defining collimator and a series of four quadrupoles designed to produce focusing of the 
beams towards the detector. Finally the two beams are again split in the vertical plane 
and recombined in a second achromat, where three stations of a Micromega-type [12j 



detector form a kaon beam spectrometer |13j (not used in the present analysis). 

After passing the cleaning and the final collimators, the beams enter the decay volume 
housed in a 114 m long cylindrical vacuum tank with a diameter of 1.92 m for the first 
65 m, and 2.4 m for the rest. Both beams follow the same path in the decay volume: 
their axes coincide within 1 mm, while the transverse size of the beams is about 1 cm. 
With 7 X 10^^ protons incident on the target per burst of ~ 4.8 s duration, the positive 
(negative) beam flux at the entrance of the decay volume is 3.8 X 10^ (2.6 X 10^) particles 
per pulse, of which 5.7% (4.9%) are The / K~ flux ratio is about 1.8. The 

fraction of beam kaons decaying in the decay volume at nominal momentum is 22%. 

The decay volume is followed by a magnetic spectrometer used for reconstructing 
TT^TT^TT decays. The spectrometer is housed in a tank filled with helium at nearly 
atmospheric pressure, separated from the vacuum tank by a thin (0.31%Xo) Kevlar com- 
posite window. A thin-walled aluminium beam pipe of 16 cm outer diameter traversing 
the centre of the spectrometer (and all the following detectors) allows the undecayed 
beam particles and the muon halo from decays of beam pions to continue their path 
in vacuum. The spectrometer consists of four drift chambers (DCH): DCHl, DCH2 lo- 
cated upstream, and DCH3, DCH4 downstream of a dipole magnet. The magnet has a 
field integral / Bydz = 0.4 Tm, thus providing a horizontal transverse momentum kick 
APj. = 120 MeV/c for charged particles. The DCHs have the shape of a regular octagon 
with a transverse size of about 2.8 m and a fiducial area of about 4.5 m^. Each chamber 
is composed of eight planes of sense wires arranged in four pairs of staggered planes ori- 
ented horizontally, vertically, and along each of the two orthogonal 45° directions. The 
spatial resolution of each DCH is ax = Cy = 90yum. The nominal momentum resolution 
of the magnetic spectrometer is parameterized as CTp/p = (1.02 © 0.044 ■ p)% {p expressed 
in GeV/c). The measured resolution on the reconstructed 37r^ invariant mass varied 
during the running period corresponding to the considered data sample in the range of 
(1.65 — 1.72) MeV/c^, depending on DCH performance. 

The magnetic spectrometer is followed by a plastic scintillator hodoscope (HOD) used 
to produce fast trigger signals and to provide precise time measurements of charged par- 
ticles. The hodoscope has a regular octagonal shape with a transverse size of about 2.4 
m. It consists of a plane of horizontal and a plane of vertical strip-shaped counters. Each 
plane consists of 64 counters arranged in four quadrants. Each quadrant is subdivided 
into 4 segments taking part in the trigger logic. Counter widths (lengths) vary from 6.5 
cm (121 cm) for central counters to 9.9 cm (60 cm) for peripheral ones. 

The hodoscope is followed by a liquid krypton electromagnetic calorimeter (LKr), a 
hadronic calorimeter (HAC) and a muon detector (MUV), all of which are not used in 
the present analysis. A detailed description of the components of the NA48 detector can 
be found elsewhere [S]. 

The K"^ 7r^7r+7r~ events are triggered with a two-level system. At the first level 
(LI), the rate of ~ 500 kHz is reduced to ~ 100 kHz by requiring coincidences of hits 
in the two planes of the HOD in at least two of the 16 non-overlapping segments. The 
second level (L2) is based on a hardware system computing coordinates of hits from DCH 
drift times, and a farm of asynchronous microprocessors performing fast reconstruction of 
tracks and running a selection algorithm, which requires at least two tracks to originate 
in the decay volume with the closest distance of approach less than 5 cm. LI triggers 
not satisfying this condition are examined further and accepted nevertheless if there is 



a reconstructed track which is not kinematically compatible with a tt^tt*^ decay of a 
having momentum of 60 GeV/c directed along the beam axis. The resulting trigger rate 
is about 10 kHz. 

2 Data analysis 

Reconstruction and selection 

Event reconstruction is based entirely on the magnetic spectrometer information. Tracks 
are reconstructed from hits in DCHs using the measured magnetic field map rescaled 
according to the recorded value of electric current in the spectrometer analyzing mag- 
net. Three-track vertices, compatible with a K"^ — »• 7r^7r"'"7r~ decay, are reconstructed by 
extrapolation of track segments from the upstream part of the spectrometer back into 
the decay volume, taking into account the stray magnetic fields due to the Earth's field 
and parasitic magnetization of the vacuum tank, and multiple scattering in the Kevlar 
window. The stray field correction is based on a three-dimensional field map measured in 
the entire vacuum tank, and reduces the amplitude of the observed sinusoidal variation 
of the reconstructed Svr invariant mass on the azimuthal orientation of the odd pion by 
more than an order of magnitude to a level below 0.05 MeV/c^. The event kinematics 
is calculated using measured momenta and track directions extrapolated to the decay 
vertex. 

The principal selection criteria applied to the reconstructed variables are the following: 

• Total charge of the three pion candidates: Q = ±1; 

• Transverse momentum: Pt < 0.3 GeV/c; 

• Longitudinal vertex position within the decay volume: Z-^tx > Z final coii.', 

• Transverse vertex radius within the beam area: R^tx < 3 cm; 

• Kaon momentum within the nominal range: 54 GeV/c < \Pk\ < 66 GeV/c. 

To improve the resolution on the kinematic variables, and to reduce the impact of 
differences between data and MC resolutions, the events were passed through a kinematic 
fitting procedure with three constraints (constraining the initial kaon direction to by along 
the z axis, and the invariant mass to the kaon mass). Events with the quality of the 
kinematic fit corresponding to probability p < 10"^ were rejected. The fraction of these 
rejected events increases as a function of deviation of the reconstructed Sir mass from the 
PDG kaon mass |AM|; in particular, all the events with |AAf | > 10 MeV/c^ are outside 
the signal region, i.e. rejected by the above condition. 

The geometric acceptance for the tt^tt^tt^ decays is mainly determined by the 

beam pipe traversing the centres of the DCHs, and the material in the central region 
of each DCH where central DCH wires terminate. This material defines a region of high 
DCH inefficiencjl^. This inefficiency together with beam optics performance and variations 
influences the acceptance, and is difficult to reproduce accurately with a MC simulation. 

^The acceptance is not biased by the finite outer size of the DCHs due to a relatively small Q- value 
of the -> 7r±7r+7r- decay: Q = 75.0 MeV. 



To minimize the effects of this problem, it is required that the transverse positions of each 
pion in DCHl and DCH4 planes {i = 1,4) satisfy the condition \R.,^i — i?o| > 18 cm, 
where is the position of the momentum-weighted average of the three pions' impact 
points: Rq = YH^iiRnilPml) / iP-Kil (for DCH4 plane, trajectories of pions are linearly 
extrapolated from DCHl and DCH2 planes). Rq corresponds to the transverse position 
of the line of flight of the initial kaon. The value of 18 cm was chosen to exclude safely 
the inefficient central region taking into account the beam sizes and variations of their 
average transverse positions. The described selection criterion costs about 50% of the 
statistics; however an appropriate MC description of the experimental conditions is more 
important than the sample size for the present analysi^. As will be shown below in the 
discussion leading to the results presented in Fig. [3l the MC simulation of the experimental 
conditions reproduces the data distributions to a level of a few parts per mille. 

The selection leaves a sample of 4.71 x 10^ events, which is practically background 
free, as — > vr^vr+Tr" is by far the dominant decay mode of the charged kaon with 
more than one charged particle in the final state. The fact that backgrounds due to other 
decays of beam kaons and pions are negligible was also checked with a MC simulation. 

The distribution of the reconstructed Stt^ invariant mass of data events (before the 
kinematic fitting) and its comparison with MC are presented in Fig. [2^. The non-Gaussian 
tails of the mass distribution are primarily due to vr^ —>■ ^^v^i decays in flight, and are 
well understood in terms of MC simulation. The ratio of MC to data mass spectra is 
presented in Fig. [2)d. It demonstrates the imperfection of resolution description in MC, 
and a deficit of MC events in the low mass region, which however contains a small fraction 
of the events, and is mostly outside the signal region. The Dalitz plot distribution of the 
selected data events (after the kinematic fitting) F^ataiu, \v\) used for the subsequent 
analysis is presented in Fig. The bin sizes of the Dalitz plot distributions used in the 
analysis are 6u = 5v = 0.05. 

Correction for trigger inefRciency 

To simplify the treatment of the trigger inefficiency, stable trigger performance was the 
main condition used to select the sample to be used for the analysi^. Inefficiencies of 
both LI and L2 trigger components were directly measured as functions of {u, \v\) using 
control data samples of prescaled low bias triggers collected along with the main triggers, 
which allowed a correction of the observed {u, \v\) distributions, and propagation of the 
statistical errors of trigger inefficiencies into the result. 

The LI trigger condition requiring a coincidence of hits in two of the 16 non-overlapping 
HOD segments is loose, as there are three charged particles in a fully reconstructed event, 
providing a rather low (and stable in time) inefficiency of 0.6 x 10~^ for the selected event 
sample. However, the LI inefficiency depends rather strongly on the kinematic variables. 
The primary mechanism generating such a dependence is the enhancement of inefficiency 
for topologies with two pions hitting the same HOD segment; such events preferably 
belong to the kinematic regions characterized by small relative velocity of a certain vrvr 
pair in the kaon rest frame. 

The L2 inefficiency, which is due to local DCH inefficiencies affecting the trigger more 

^In a different case, the charge asymmetry analysis [5] is performed with soft cuts maximizing the 
selected data sample, and is based on cancellation, rather than simulation, of the systematic effects. 
^As it was already noted, the size of the data sample is not a limitation for this analysis. 
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Figure 2: (a) Reconstructed spectrum of Svr^ invariant mass M(37r) (upper line) and its 
presentation in terms of normalized MC components: (1) events without it —>■ fiu decay 
in flight, (2) events with it ^ fiu decay, (3) IB radiative K^t^^ events, (b) The ratio of 
MC/data M(37r) spectra demonstrating the imperfection of resolution description in MC, 
and a deficit of MC events at low M(37r) outside the signal region, (c) Reconstructed 
distribution in the kinematic variables F^ataiu, \v\) (after kinematic fit). 



strongly than the off-line reconstruction due to lower redundancy and trigger timing 
effects, was measured to be (0.32 ± 0.05) x 10^^ (the error indicates the maximum size of 
its variation during the data taking period). It did not exhibit any significant correlation 
to the kinematic variables, due to relative complexity of the decision-making algorithm. 

Monte Carlo simulation 

A detailed GEANT-based MC simulation was developed, which includes full detector 
geometry and material description, simulation of stray magnetic fields, DCH local in- 
efficiencies and misalignment, the beam line (which allows a reproduction of the kaon 
momentum spectra and beam profiles), and / K~ relative fluxes. Moreover, time vari- 
ations of the above effects during the running period were simulated. A production of 
6.7 X 10^ K"^ —>■ TT^TT+TT^ events distributed according to the matrix element ([2]) with PDG 
values of the slope parameters [Ij was performed to determine the detector acceptance. 
A sample of 1.16 x 10^ MC events (almost 2.5 times larger than the data sample) passes 
the selection. 

Comparison of data and MC distributions in such significant variables as longitudinal 
decay vertex position and illuminations of DCHl and DCH4 planes by pions is presented 
in Fig. [31 and demonstrates that MC simulation reproduces the data distributions to a 
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Figure 3: Left column: reconstructed data distributions, right column: the corresponding 
ratios of data/MC distributions of (a) vertex z position, and pion radial position (for each 
of the 3 pions) in the planes of (b) DCHl and (c) DCH4. 



level of a few units of 10 ^. The precision of the data description can be improved by 
tighter cuts on pion radial positions. 

Fitting procedure 

The measurement method is based on fitting the binned reconstructed data distribution 
Fdataiu, \v\) presented in Fig. [2b with a sum of four reconstructed MC components gener- 
ated according to the four terms in the polynomial ([2]) presented in Fig. HI Let us denote 
these reconstructed MC distributions as -Fo(m, |f |), |f |), FuZ^u, \v\), and Fy2(u, \v\). 

To obtain them, the MC sample (distributed in kinematic variables according to the PDG 
slope parameters) was divided into four subsample£], and events in each subsample were 
assigned appropriate weights depending on the generated (m, |f |) to obtain the desired 
distributions. 

The following functional defining the agreement of the shapes of data and MC distri- 

■^The relative sizes of the four subsamples were subject of optimization in order to minimize the 
statistical error of the measurement. The F^^iu, \v\), i^„2(u, \v\), and _F!„2(u, samples are of equal sizes, 
while the Fq{u, \v\) sample is 5 times larger than each of the former three. 




Figure 4: Reconstructed distributions in the kinematic variables (M,|f |) of the four MC 
components: (a) Fo{u, \v\), (b) \v\), (c) F„2(m, |f |), (d) F^2{u, \v\). 



butions is minimized using the MINUIT package [15] in order to measure the values of 
the slope parameters {g, h, k): 



The sum is evaluated over all the bins of reconstructed {u, \v\) distributions with at least 
1000 data events, which eliminates the necessity to include non-Gaussian behaviour of 
errors. Here 6Fdata{u, \v\) is an uncertainty of the number of data events in a given 
bin (composed of a statistical part and a trigger efficiency part added in quadrature), 
FMc{g, h, k, u, \v\) is a MC population of a bin for given values of {g, h, k): 

FMc(,g,h,k,u,\v\) = Fo{u,\v\)/Io + gFu{u,\v\)/Iu+ (6) 

hFu2{u, \v\)/Iu^ + kF^2{u, \v\)/Iy2, 

and 6Fmc{9, h, k, u, \v\) is its statistical error: 

6'FMci9,h,k,u,\v\) = 6'Foiu,\v\)/I^ + g'6'F^iu,\v\)/ll+ (7) 



Parameter 


Value 


Uncertainties 


Trigger correction 


statistical MC stat. 


g X 10^ 
h X 10^ 
k X 10^ 


-21.134 ± 0.013 
1.848 ± 0.022 
-0.463 ± 0.007 


0.009 0.008 
0.015 0.013 
0.005 0.004 


-0.008 ± 0.005 
0.116 ± 0.009 
0.033 ± 0.003 



Table 1: Fit results with statistical and MC statistical uncertainties, trigger corrections 
and their statistical uncertainties. The trigger corrections are included into the resulting 
values. 

h''6^F^2{u, \v\)lll + k^d^F^u, \v\)lll,. 

Here /q, /«, and 7^2 are the normalization constants computed taking into account the 
numbers of generated events in each of the four MC subsamples, and the integrals of the 
four terms in ([2]) over the Dalitz plot. The free parameters of the functional (jSD are the 
slope parameters [g, h, k) and an overall MC normalization parameter A^. 

The minimization yields x^/NDF = 1669/1585, corresponding to a satisfactory prob- 
ability of 7.0%. The results of the fit and the trigger corrections are presented in Tabled! 
The non-zero values of the corrections arise mostly from the LI trigger inefficiency depen- 
dence on kinematic variables, while their statistical errors receive contributions of similar 
sizes from LI and L2 trigger efficiency uncertainties. 

Keeping only the linear term gu in ([2]) yields a fit of unacceptable quality: x^/NDF = 
13683/1583. Including the terms proportional to and uv"^ (the only cubic terms allowed 
by Bose symmetry) yields values for the corresponding cubic slope parameters compatible 
with zero. 

Stability checks 

Stability of the results with respect to variations of the selection conditions on vertex fit 
quality, Pt, Rvtx, \Pk\ and R^, and binning variations was checked. Stability with respect 
to exclusion of (m, |f |) bins with large deviations of the Coulomb factor ([3]) from unity and 
with respect to kaon sigij^ was checked as well. No statistically significant dependencies 
were found. Stability with respect to various ways of binning the data was checked; 
comparison of slope measurements in bins of reconstructed longitudinal coordinate of 
decay vertex Z^tx (since the acceptance depends strongly on this variable), and in data 
taking periods are shown in Fig. [5] as the most significant examples. 

Systematic uncertainties 

The Coulomb factor ([3]) used in description of the event density contains a pole C{u, v) — >• 
oo corresponding to a pair of opposite sign pions having zero relative velocity. The 
implemented Monte Carlo simulation involves a certain approximation to C{u,v) in the 
pole region. In view of that, sensitivity of the result to treatment of the pole was studied, in 
particular, by using an alternative fitting method involving projections of the distributions 
in each of kinematic variables. The assigned conservative systematic uncertainties due to 
the fitting procedure are listed in Table [2j They are expected to be diminished to a 
negligible level in a future analysis of the full statistics. 

^Combined and K~ sample is used to obtain the result. Stability of the slope parameters with 
respect to kaon sign is a consequence of the experimental fact [5J that CP invariance holds at the discussed 
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Figure 5: Slope measurements in bins of Z^tx (left column) and in running periods (right 
column). Systematic errors are not shown. 



An important source of systematic uncertainty is the imperfect description of the 
resolution in pion momentum, which can be observed in Fig. [2]d as a slight disagreement 
of the shapes of the reconstructed M(37r) spectra for data and MC in the signal regioij^. 
To evaluate the corresponding effect, two different plausible ways of introducing smearing 
of the MC resolution were used: either by increasing the smearing of DCH space points 
from 90/im to 100/im, or by adding an extra 0.09%Xo layer of matter in the position 
of the Kevlar window (the former correction is more realistic). The sizes of the added 
perturbations are such as to correct for the "double bump" shape of the ratio in the signal 
region. Both methods lead to similar systematic uncertainties on the slope parameters 
attributed to description of resolution in pion momentum. These uncertainties are listed 
in Table H 

Effects due to imperfect knowledge of the magnetic field in the spectrometer magnet 
were evaluated. The variation of the magnet current can be monitored with a relative 
precision of 5 x 10"*^. Smaller variations are continuously controlled with a precision of 



level of precision. 

^As was discussed above, a large disagreement at low M(37r) values is outside the signal region. 



Effect 


5g X 10^ 


5h X 10^ 


5k X 10^ 


Fitting procedure 


0.009 


0.007 


0.006 


Pion momentum resolution 


0.004 


0.031 


0.009 


Spectrometer magnetic field 


0.002 


0.008 


0.004 


Spectrometer misalignment 


0.002 


0.002 


0.001 


Stray magnetic field 


0.001 


0.002 


0.001 


Total 


0.010 


0.033 


0.012 



Table 2: A summary of the systematic uncertainties. 



~ 10~^ by the deviation of the measured charge-averaged kaon mass from the nominal 
PDG value. A time-dependent correction is introduced by scaling the reconstructed pion 
momenta, decreasing the effect of overall field scale to a negligible level. To account for 
possible differences between the shape of the field map used for simulation and the true 
field, variations of the MC field map were artificially introduced, consistently with the 
known precision of field measurement of ~ 10~^. The corresponding uncertainties are 
listed in Table [21 

The transverse positions of DCHs and individual wires were controlled and realigned 
at the level of reconstruction software every 2-4 weeks of data taking using data collected 
in special runs in which muon tracks were recorded with no magnetic field in the spec- 
trometer. This allows an alignment precision of ~ 30 /xm to be reached. However, time 
variations of DCH alignment on a shorter time scale can bias the measurement. These 
variations were measured by the difference between the average reconstructed 37r invariant 
masses for and K~ decays, and taken into account. The precision with which these 
effects are simulated leads to systematic uncertainties presented in Table HI 

The effects due to limited precision of measurement of the stray magnetic field in the 
decay volume were estimated by variation of the stray field map used for decay vertex 
reconstruction; the corresponding systematic effects are presented in Table |2j 

The kaon momentum spectra were carefully simulated, and the related residual uncer- 
tainties were found to be negligible. Possible differences between data and MC transverse 
scales were found to have a negligible influence on the result. 

The total systematic errors were obtained by summing the above contributions in 
quadrature, and are presented in Table [2l 

Conclusions 

The Dalitz plot slope parameters of the Tr^Tr+Tr" decays measured with a fraction 

of NA48/2 data sample ignoring radiative effects (apart from the Coulomb factor) and 
strong rescattering effects, are: 

^ = (-21.134 ±0.017)%, /i = (1.848 ± 0.040)%, = (-0.463 ± 0.014)%. 

These values are in agreement with the world average^, and have an order of magnitude 
smaller uncertainties. This is the first measurement of a non-zero value of the quadratic 

'^The PDG averages results separately for and K~ [T]; averaging the PDG data between and 
K~ decays should take into account correlated systematic uncertainties of the and measurements 
by the same experiment [2J. 



slope parameter h. The compatibility of the measured distribution with the PDG polyno- 
mial parameterization [I] appears still to be acceptable at an improved level of precision; 
no significant higher order slope parameters were found. 

The measurement of the slope parameters is in agreement with a full next-to-leading 
order computation [T6] : 



^ = (-22.0 ±2.0)%, /i = (1.2 ±0.5)%, A; = (-0.54 ±0.15)%. 

The whole NA48/2 sample suitable for K"^ Svr^ Dalitz plot distribution analysis 
contains at least three times more data; a more elaborate analysis is foreseen when the 
corresponding theoretical framework is available. 
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